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ABSTRACT 

We present the catalog of X-ray sources detected in a shallow Chandra survey of the inner 2°x0.8° 
of the Galaxy, and in two deeper observations of the Radio Arches and Sgr B2. The catalog contains 
1352 objects that are highly-absorbed (iVH>4 x 10^^ cm~^) and are therefore likely to lie near the 
Galactic center (iD«8 kpc), and 549 less-absorbed sources that lie within <6 kc of Earth. Based 
on the inferred luminosities of the X-ray sources and the expected numbers of various classes of 
objects, we suggest that the sources with Lx ^ 10-^^ erg s~^ that comprise ~90% of the catalog 
are cataclysmic variables, and that the «100 brighter objects are accreting neutron stars and black 
holes, young isolated pulsars, and Wolf-Rayet and O stars in coUiding-wind binaries. We find that 
the spatial distribution of X-ray sources matches that of the old stellar population observed in the 
infrared, which supports our suggestion that most of the X-ray sources are old cataclysmic variables. 
However, we find that there is an apparent excess of ~10 bright sources in the Radio Arches region. 
That region is already known to be the site of recent star formation, so we suggest that the bright 
sources in this region are young high-mass X-ray binaries, pulsars, or WR/0 star binaries. We briefly 
discuss some astrophysical questions that this catalog can be used to address. 
Subject headings: catalogs — Galaxy: center — X-rays: general 



I 1. INTRODUCTION 

The exquisite sensitivity of the Chandra X-ray Ohser- 

' vatory has provided vast improvements in our under- 
standing of faint, hard X-ray sources. Chandra obser- 
vations of distant galaxies allow us to study the X-ray 

I population at luminosities similar to those accessible in 
our own Galaxy with wide-field X-ray instruments like 
the Rossi X-ray Timing Explorer All-Sky Monitor and 
the BeppoS AX Wide-Field Camera (Lx ~ 10^^ erg s'^), 
while observations of our own Galaxy are sensitive to 
sources a million times fainter than found with previ- 

I ous wide-field surveys. One of the most dramatic prod- 
ucts of this improvement in sensitivity are the Chan- 
dra observations of the Galactic center (Wang, Gotthelf, 

, & Lang 2002; Baganoff et al. 2003). Whereas previ- 
ous imaging surveys identified dozens of X-ray sources 
against a back g round of bright Gala ctic diffuse emission 
llWatson et alJ I198H fPavhnskv. G rebenev fc Sunva^ 
1994'"Prcdchl & Trucmocr 1995; Sid oh et alJl200iLfl9M 
Sakano et al. 2002), Chandra observations h ave revealed 
thousands of individual X-ray sources fe.g.. lWang et al.l 
l2002t iMuno et al.ll200 j|) and d iscret e, fi lamenta ry fea- 
tures (e.g., iLu. Wang, fc Lang|l2003t mTscf-Zadc h et all 
12005a'). This high concentration of X-ray sources is not 
surprising. The 2°x0.°8 (300x125 pc for a distance of 
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8 kpc; iMcNamara et al .' '2000) region centered on the 
Galactic center contains roughly 1% of the Galactic mass 
(Lauhardt, Zylka, & Mezger 2002). Moreover, unlike the 
Galactic bulge, star formation has occurred continuously 
in the central region of the Galaxy, as is strikingly il- 
lustrated by the >60 ultra-co mpact HH regions in the 
giant molecular cloud Sgr B2 l|de Pree. Goss. fc Gaumd 
Il998j) . and by three young, dense clusters of massive stars 
(the Arches, the Quint uplet, and t he Central Parsec; 
iKrabbe et al] IT995: Fige r et al.l [l999'). 

A wealth of questions can be addressed with Chan- 
dra observations, and with subsequent comparisons to 
multi-wavelength catalogs. For instance, a variety of 
studies of the synthesis of compact, accreting binaries 
have been designed to explain the large number of X- 
ray sources Chandra detects in the Galactic center, and 
the results constrain, for example, the amount of angu- 
lar momentum dissipated in the common envelope phase 
(e.g., Pfahl, Rappaport, & Podsiadlowski 2002; Belczyn- 
ski & Taam 2004; Liu & Li 2005; Ruiter, Belczynski, fc 
Harrison 2005). Chandra also can detect outbursts from 
transient low-mass X-ray binaries (LMXBs) at much 
lower flux levels than are accessible with traditional wide- 
field X-ray surveys, which provides unique insight into 
the duty cycles and emission mechanisms of compact 
objects accreting at v ery low rates (M < lO"^ " !" Mp^ 
vr-^: e.g..lKing,2000l^^yiinands■ Miller, fc Wand 12002 
Wiinands fc Wand I2002t lin't Zandl 120051 ISakano eTall 
20051) . Combining Chandra and radio observations 
can be used to identify young stars with powerful 
winds, which helps to constrain the rate at which mas- 
sive stars have formed recently in the Galactic cen- 
ter HYusef-Zadeh et al.ll2f)0l ILaw fc Yusef-Zadehll200l 
IMuno et alJl2005bD . 

Several observational studies have discussed the popu- 
lation o f X-ray s ources in Chandra surveys of the Galactic 
center. I Wang et al. (2002.) presented images from shal- 
low (12 ks) Chandra exposures of the 2°x0.°8 around 
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Fig. 1. — Mosaic image of the 0.°8x2° field along the Galactic plane, centered on Sgr A*. The raw image has been adaptively smoothed 
with the CIAO tool csmooth for display purposes. The prominent features are the Sgr A complex at the center of the image, two bright 
X-ray binaries and their dust scattering halos at (=0.275°and —0.88°. Only a fraction of the brightest point sources are visible in this 
image. Note that the smoothing algorithm introduces significant artifacts, especially at the edges of the deep observations centered on Sgr 
B2 and the Arches. 



Sgr A*, and gave a general overview of the number 
of X-ray sources and the properties of the diffuse X- 
ray emission. Takagi, Murakami, & Koyama (2002) 
studied the properties of X-ray sources asso ciated with 
HII regions in Sgr B2. lYuseif-Zadeh et all (|2002) and 
iLaw fc Yusef-Zadebl l)2004j) and studied the X -ray emis- 
sion from massive stars in several clusters, including 
the Arches and Quintuplet. Finally, Muno et al. (2003, 
l20'04bl l2005aj) . presented a comprehensive study of the 
population of X-ray sources with Lx = lO'^^ — 10^"^ erg 
s^^ that were discovered in a deep (625 ks) set of Chandra 
observations of the inner 25 pc around Sgr A*. However, 
aside from the inner 25 pc, no catalog containing all of 
the X-ray sources found within the inner 150 pc of the 
Galactic center has been published. 

In this paper, we rectify this by reporting the loca- 
tions and basic properties of the X-ray sources detected 
in Chandra observations of the inner 2°x0.°8 around the 
Galactic center. In >i2.1\ we present the locations of these 
X-ray sources (exc luding the inner 8' covered by the cata- 
log in lMuno et alJl2b 031 in order to facilitate searches for 
multi- wavelength counterparts. In H2.2\ we present the 
fluxes and basic spectral properties in order to constrain 
the origin of the X-ray emission, and to serve as baseline 
measurements for future searches for transient sources. 
In ^'6.1\ we examine the spatial distribution of the X-ray 
sources to determine how they are related to the stellar 
population that is observed in the infrared. In ^'6.2\ we 
study the luminosity distribution of the X-ray sources, 
and report variations in the relative numbers of bright 
sources that could be related to recent star formation. 

2. OBSERVATIONS 

The central 300x125 pc of the Galaxy has been ob- 
served on several occasions with the imaging array of the 
Chandra Advanced CCD Imaging Spectrometer (ACIS-I 
IWeisskopf et al...,2002') . The entire r egion was surveye d 
with overlapping « 12 ks observations ijWang et all2002|) . 



and one additional short observation was obtained cen- 
tered on the LMXB IE 1740.7-2942. Deeper observa- 
tions were taken of the H II regions in the giant molecular 
cloud Sgr B2 (100 ks; Takagi et al."2002^, and the non- 
thermal radio features referr ed to as t he Arche s (50 ks o f 
pubhc data as of 2005 June: ILaw irY uscf-Za dehll2004j) . 
We list each of these observations in Tabled and present 
a mosaic image of the survey in Figure ^ 

The ACIS-I is a set of four 1024-by-1024 pixel CCDs, 
covering a field of view of 17' by 17'. When placed on- 
axis at the focal plane of the grazing-incidence X-ray 
mirrors, the imaging resolution is determined primarily 
by the pixel size of the CCDs, 0'.'492. The CCDs also 
measure the energies of incident photons within a cal- 
ibrated energy band of 0.5-8 keV, with a resolution of 
50-300 eV (depending on photon energy and distance 
from the read-out node). The CCD frames are read out 
every 3.2 s, which provides the nominal time resolution 
of the data. 

We reduced the observations using standard tools from 
the CIAO package, version 2.3.01. We started with the 
level- 1 event lists provided by the Chandra X-ray Center 
(CXC), and removed the pixel randomization applied by 
the default processing software. We then modified the 
pulse heights of each event to partially correct for the 
position-dependent charge-transfer inefficiency caused by 
radiation damage early in the mission, using software 
provided by Townslev et al. (20023)- We excluded most 
events flagged as possible background, but left in possi- 
ble cosmic ray afterglows because in the version of the 
processing software that we used they were difficult to 
distinguish from genuine X-rays from the strong diffuse 
emission and numerous point sources in the field. We 
applied the standard ASCA grade filters to the events, 
as well as the good-time filters supplied by the CXC. 
Finally, we searched each observation for time intervals 
when the detector background flared to > 3cr above the 
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TABLE 1 

Observations of the Central 2°x0.8°of the Galaxy 











Aim Point 




Start Time 


Sequence 


Target 


Exposure 


RA DEC 


Roll 


(UT) 






(ks) 


(degrees J2000) 


(degrees) 



2000 


Aug 30 16:59:32 


658 


IE 1740.7-2942 


9.2 


265.97583 


-29.75008 


270.8 


2000 


Mar 29 09:44:36 


944 


SGR B2 


97.5 


266.78034 


-28.44169 


87.8 


2000 


Jul 


07 


19:05:19 


945 


GC ARC 


48.8 


266.58192 


—28.87196 


284.4 


2001 


Jul 


19 


10:01:48 


2267 


GCS 


20 


8.7 


266.17150 


-29.27337 


283.8 


2001 


Jul 


20 


04:37:11 


2268 


GCS 


21 


10.8 


265.98136 


—29.17141 


283.8 


2001 


Jul 


16 


02:15:50 


2269 


GCS 


1 


10.5 


267.05495 


-28.37576 


283.8 


2001 


Jul 


20 


08:00:49 


2270 


GCS 


22 


10.6 


266.24512 


-29.54138 


283.8 


2001 


Jul 


16 


05:35:55 


2271 


GCS 


2 


10.4 


266.86502 


—28.27455 


283.8 


2001 


Jul 


20 


11:12:40 


2272 


GCS 


23 


11.6 


266.05423 


—29.43957 


283.8 


2001 


Jul 


18 


00:48:28 


2273 


GCS 


10 


11.2 


266.70988 


—28.87565 


283.8 


2001 


Jul 


16 


08:44:25 


2274 


GCS 


3 


10.4 


266.67662 


-28.17301 


283.8 


2001 


Jul 


20 


14:41:10 


2275 


GCS 


24 


11.6 


265.86371 


-29.33729 


283.8 


2001 


Jul 


18 


04:16:58 


2276 


GCS 


11 


11.6 


266.51970 


-28.77438 


283.8 


2001 


Jul 


16 


11:52:55 


2277 


GCS 


4 


10.4 


266.94061 


—28.54231 


283.8 


2001 


Jul 


20 


18:09:40 


2278 


GCS 


25 


11.6 


266.12769 


-29.70775 


283.8 


2001 


Jul 


18 


07:45:28 


2279 


GCS 


12 


11.6 


266.33020 


-28.67281 


283.8 


2001 


Jul 


16 


15:01:25 


2280 


GCS 


5 


10.4 


266.75037 


-28.44124 


283.8 


2001 


Jul 


20 


21:38:10 


2281 


GCS 


26 


11.6 


265.93652 


—29.60557 


283.8 


2001 


Jul 


18 


11:13:58 


2282 


GCS 


13 


10.6 


266.59425 


—29.04216 


283.8 


9001 


Jul 


21 


01:06:39 


2283 


GCS 


27 


1 1 




on c;o'^1 ^ 


zoo.o 


2001 


Jul 


18 


14:25:48 


2284 


GCS 


14 


10.6 


266.40487 


-28.94088 


283.8 


2001 


Jul 


16 


18:09:55 


2285 


GCS 


6 


10.4 


266.56112 


-28.34029 


283.4 


2001 


Jul 


21 


04:35:09 


2286 


GCS 


28 


11.6 


266.00997 


-29.87372 


283.8 


2001 


Jul 


18 


17:37:38 


2287 


GCS 


15 


10.6 


266.21439 


-28.83925 


283.8 


2001 


Jul 


17 


14:11:51 


2288 


GCS 


7 


11.1 


266.82518 


-28.70891 


283.8 


2001 


Jul 


21 


08:03:39 


2289 


GCS 


29 


11.6 


265.81855 


-29.77165 


283.8 


2001 


Jul 


21 


11:32:10 


2290 


GCS 


30 


11.6 


265.62772 


-29.66900 


283.8 


2001 


Jul 


18 


20:49:28 


2291 


GCS 


16 


10.6 


266.47839 


-29.20880 


283.8 


2001 


Jul 


17 


17:51:28 


2292 


GCS 


8 


11.6 


266.63516 


-28.60795 


283.8 


2001 


Jul 


19 


00:01:18 


2293 


GCS 


17 


11.1 


266.28794 


-29.10740 


283.8 


2001 


Jul 


17 


21:19:58 


2294 


GCS 


9 


11.6 


266.44581 


-28.50671 


283.8 


2001 


Jul 


19 


03:21:28 


2295 


GCS 


18 


11.1 


266.09836 


-29.00518 


283.8 


2001 


Jul 


19 


06:41:38 


2296 


GCS 


19 


11.1 


266.36205 


-29.37522 


283.8 



mean level, and removed such intervals when they oc- 
curred (in ObsIDs 2267, 2269, 2273, 2288, and 944). 

2.1. Source Detection and Initial Localization 

We searched for X-ray sources separately in sets of 9 
images for each observation usi ng the wavelet routine 
wavdetect ("Freem an et al.ll2002D . We generated images 
in three energy bands: the full 0.5-8.0 keV band, the 0.5- 
2.0 keV band to increase our sensitivity to foreground 
sources, and 4-8 keV to increase our sensitivity to highly 
absorbed sources. For the purposes of source detection 
only, we removed events that had been flagged as possible 
cosmic ray afterglows. We employed the default "Mex- 
ican Hat" wavelet, and used a sensitivity threshold of 
10~^ that corresponds to chance of detecting a spurious 
source per pixel if the local background is spatially uni- 
form. We searched each energy band using a succession 
of three images centered on the aim point of each expo- 
sure: 1024x1024 images at the full Chandra resolution 
of 0'.'5, 1024x1024 images binned by a factor of 2 to a 
resolution of 1", and images of variable size that covered 
the entire ACIS-I exposure for each observation with a 
resolution of 2". We used wavelet scales that increased 
by a factor of 1-4 for the 0'.'5 image, 1-8 for the 
1" image, and 1-16 for the 2" image. This succession of 
three images and spatial scales were chosen because the 



point-spread function (PSF) broadens as a function of 
offset from the aim point. The resulting source list from 
the 33 observations contained 1901 unique sources, 225 
of which were only detected in the soft band, and 382 of 
which were only detected in the hard band. Based on the 
sensitivity threshold for wavdetect (10~^), we expect 2 
spurious sources per field, or ~70 in the entire survey. 

We attempted to refine the astromctry for each obser- 
vation by matching foreground X-ray sources detected 
in the soft ban d, many of which are likely to be K and 
M dwarf stars l)Muno et alJl2003|) . to infrared sources in 
the Two-Micron All-Sky Survey (2MASS) catalog. For 
the two observations longer than 50 ks (Sgr B2 and the 
Arches), we found that there were «20 matches between 
the soft X-ray and 2MASS catalog within 5' of the aim 
point. By randomly shifting the relative positions of the 
two catalogs, we determined that, with 90% confidence, 
fewer than 25% of these matches should be random. 
Based on these matches, we could derive the absolute 
astrometry of the Chandra pointing to within O'.'l. 

Unfortunately, the shorter observations were less sen- 
sitive, and contained far fewer X-ray sources. Whereas 
the > 50 ks exposures contained «200 X-ray sources, 
the shorter exposures generally contained only a cou- 
ple dozen. In general, <2 X-ray sources could be iden- 
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tified with 2MASS counterparts, which was insufhcient 
to improve the pipeUne astrometry. Therefore, starting 
with the fields adja cent to the deep exposures of Sgr A* 
l)Muno et all 120031) . the Arches, and Sgr B2 and mov- 
ing outward, we derived the astrometry by matching the 
wavdetect positions of the X-ray sources matched those 
of adjacent fields. Unfortunately, because most of the 
matches made in this way relied on X-ray sources de- 
tected > 6' off-axis, the statistical uncertainty on the 
wavdetect positions of each X-ray sources was signifi- 
cant, and the corrected astrometry again was not signif- 
icantly more accurate than the default pipeline values. 
Therefore, we expect that the astrometry will on aver- 
age be accurate to 0'.'7, but that «10% of the pointings 
will have systematic errors of up to 1". Unfortunately, 
there is no way of knowing a priori which observations 
have the larger astrometric errors, so our astrometric ac- 
curacy is limited to «1" in the shallower exposures. 

In Table IBll we list the refined positions and posi- 
tional uncertainties (90% confidence) , which include sta- 
tistical and systematic terms. For the statistical posi- 
ti onal error, we use the c ount-rate-dependent estimates 
in [Alexander et all 1)200,'^ . For the majority of sources 
with <200 net counts, we assume the uncertainty A=0'.'6 
for offsets d<5' and 

A = 0.6+ (^^^) arcsec (1) 

for 9>5\ For 19 bright sources with >200 net counts, we 
assume A=0'.'3 for offsets 9<5' and 

A = 0.3 + (-^^^ arcsec (2) 

for 6>5\ We combine the statistical errors in quadra- 
ture with 0'.'7 systematic uncertainties for the shallow 
observations (all those except ObsIDs 944 and 945, for 
which we assume no systematic uncertainty). In Ta- 
ble IBll we also list the observation in which a source 
was detected in, the offset from the aim point of that 
observation, and the total live time each source was ob- 
served with. We omitted two bright, previously-known 
LMXBs from the table (IE 1740.7-2942 and 2E 1743.1- 
2842), because they were badly saturated in our survey 
images, and their positions and properties are reported 
elsewhere fe-g.- IMarti et al.ll2000l: IPorauet et al.ll2003(l . 

2.2. Photometry 

We computed photometry for each source in the 0.5- 
8.0 keV band using the acis_extract routine from the 
Tools for X-ray Analysis (TARA).^ We extracted event 
lists for each source for each observation, using an ex- 
traction region designed to enclose a large fraction of the 
PSF. We used a PSF at a fiducial energy of 1.5 keV for 
sources detected only in the soft band (flagged with a 
'f to indicate a possible foreground source in Tab. IBl() . 
while we used a larger extraction area corresponding to 
a PSF for 4.5 keV photons when sources were detected 
in the full or hard bands. 

In most cases, we chose polygonal regions that matched 
the contours of 90% encircled energy from the PSF. How- 
ever, if the 90% contours of the PSFs of two nearby 

^ www . astro .psu. edu/xray/docs/TARA/ 



sources overlapped, we generally used a region that corre- 
sponded to a smaller fraction of the PS F. W e have flagged 
these sources as confused ('c') in Table IbTI The smallest 
extraction region that we used matched the 70% encir- 
cled energy contour. However, because the PSF grows 
significantly as a function of off-axis angle, in many cases 
a source that appeared isolated when it was located on- 
axis in one observation was indistinguishable from its 
neighbors in images from adjacent, overlapping point- 
ings. We found that no reasonable fraction of the PSF 
would isolate counts from sources more than 7' off-axis 
whose 90% encircled-energy radii overlapped those of 
their neighbors, so we did not extract photometry from 
observations in which confused sources were >7' off-axis. 
Fortunately, these sources always lay near the aim point 
of another observation (or else they could not have been 
detected in the first place), so we were still able to obtain 
photometry for these sources. 

For each source and each observation, a background 
event list was extracted from a circular region centered 
on the point source, excluding from the event list counts 
in circles circumscribing the «92% contour of the PSF 
around any point sources. We found that this value 
struck a good balance between excluding counts from 
point sources and leaving enough counts in the image to 
determine the background. The sizes of the background 
regions were chosen such that they contained approxi- 
mately 100 events for each observation. Less than 1% 
of the counts in the background regions were from point 
sources. We also computed the effective area function 
(ARF) at the position of each source for each observa- 
tion using the CI AO tool mkarf . 

The source and background event lists were used to 
compute photometry for each source in five energy bands: 
0.5-8.0 keV (the fuU band), 0.5-2.0 keV, 2.0-3.3 keV, 
3.3-4.7 keV, and 4.7 -8.0 keV. Th e se ba nds are identi- 
cal to those used by iMuno et alJ 12P03') for the cata- 
log of X-ray sources within 25 pc of Sgr A*. The en- 
ergy bands were chosen so that they sampled regions of 
the ARF with roughly constant areas, and so that the 
three high energy bands each contained about one-third 
of the net counts from most sources. The net counts in 
each band were computed by subtracting the estimated 
background from total counts. The 90% uncertainty in 
the net counts in each of the five bands were computed 
through a Bayesian analysis, with the simplifying as- 
sumption that the uncertainty on th e background was 
negligible (jKraft. Burrows, fc Nousek.1991) . When the 
90% confidence interval on the net counts was consis- 
tent with 0, we considered the upper bound on the 90% 
confidence interval to be the upper limit. 

A histogram of the number of sources as a function 
of net counts in the full band is displayed in Figure EJ 
Most sources were detected with 10-20 counts in two 
overlapping 12 ks observations. Therefore, most sources 
were only observed with «5 counts in each of the smaller 
energy bands. The net counts from 303 sources are con- 
sistent with zero at the 90% level. About 70% of these 
sources are flagged in Table IBll because their photom- 
etry could be unreliable for various reasons: some are 
detected only in the soft band where the background is 
lower, some are variable so that the mean flux is not 
meaningful, and some are confused with nearby sources. 
The remaining ss90 sources are probably spurious, as we 
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Fig. 2. — Top panel: The distribution in net counts from individ- 
ual sources. No correction is applied to account for the exposure 
across the survey, which varies by a factor of 5. Bottom panel: 
The distribution of fluxes from individual sources, derived by di- 
viding the net count rates by the effective area and exposure in 
four energy bands (0.5-2.0 keV, 2.0-3.3 keV, 3.3-4.7 keV, and 4.7- 
8.0 keV), and summing the result. There are two peaks, because 
the deeper observations were more sensitive to faint sources. In 
both panels, the solid line is used for sources located near or be- 
yond the Galactic center (iTi?>— 0.175), and the dashed line for 
foreground sources (iTi?<— 0.175). The arrows denote the median 
sensitvity for the shallow survey (labeled WGL02) and for the deep 
Arches and Sgr B2 fields. 



expected «70 spurious sources based on our detection 
threshold. We suspect we detected ~20 more spurious 
than we expected because the background in our obser- 
vations is «10 times larger than in the observations taken 
at high Galactic latitude with which wavdetect was cali- 
brated. We keep these potentially spurious sources in the 
catalog for completeness. The net counts and 90% uncer- 
tainties (or upper limits) are listed in Table IbTI These 
values are used to compute the probability of detecting 
each source given its location and exposure, which is also 
listed in the table. The Monte Carlo simulations used to 
estimate this probability are described in Appendix A. 

We computed approximate photon fluxes (in units of 
ph cm~^ s~^) for each source by dividing the net counts 
in each sub-band by the total live time (units of s) and 
the mean value of the ARF in that energy range (units of 
cm^; note that this value incorporates variations in ex- 
posure due to chip gaps and dead columns). The photon 
fluxes in the 0.5-8.0 keV energy band used throughout 
the paper are the sums of those in the sub-bands, using 
negative values when they occur (not the upper limits). 
They are listed in Table IbTI We have found that the ap- 
proximate photon fluxes that we computed differed from 
those derived from spectral flts by little more than the 



Poisson uncertainty in the count rate, because the en- 
ergy bands sampled th e ARF for the ACIS-I detector 
weU ijMuno et alJl2003|) . 

A histogram of the number of sources as a function 
of the 0.5-8.0 keV photon flux is presented in the bot- 
tom panel of Figure [3 Galactic center sources are in- 
dicated with the solid line, and foreground sources with 
the dashed line. Sources are detected with average pho- 
ton fluxes as low as 2 x 10~® ph cm~^ s~^. The largest 
number of Galactic center sources are detected near 
2 X 10"^ ph cm~2 (2.0-8.0 keV), and the largest 
number of foreground sources are found near 1 x 10~^ ph 
cm-2 s-i (0.5-2.0 keV). 

We used the counts in each energy band to compute 
two hardness ratios, which we used to characterize the 
absorption column toward each source and the steepness 
of the high-energy portion of each spectrum. The ratios 
are defined as the fractional difference between the count 
rates in two energy bands, {h — s)/{h + s), where h and 
s are the numbers of counts in the higher and lower en- 
ergy bands, respectively. The resulting ratio is bounded 
by —1 and +1. The soft color is defined by the frac- 
tional difference between counts with energies between 
2.0-3.3 keV and 0.5-2.0 keV, and the hard color using 
counts between 4.7-8.0 keV and 3.3-4.7 keV. The hard- 
ness ratios are listed in Table IbTI with uncertainties cal- 
culated according to Equation 1.31 in Lyons (1991; page 
26). 

We use the soft color {HRO) to identify foreground 
and highly-absorbed sources. By comparing the spectral 
fits and hard colors in ( Muno et a l. 2003, 2004bJ, we find 
that sources with HRO > —0.175 or that are not detected 
below 3.3 keV have absorption columns A^h > 4 x 10^^ 
cm^^, and are therefore likely to lie at or beyond the 
Galactic center. We refer to these as the "Galactic center 
sources," of which there are 1350. Sources with HRO < 
—0.5 have A^h ^ 10^^ cm~^ and he within 2 kpc of Earth. 
Sources with intermediate soft colors lie between 2-6 kpc 
from the Earth. The 549 sources with HRO < -0.175 
that are likely to be closer than 6 kpc are considered 
foreground sources. 

In Figure 121 we plot the hard color versus the flux from 
each source. Foreground sources are indicated with open 
circles, and sources at or beyond the Galactic center with 
filled circles. There are 785 Galactic center sources and 
39 foreground sources with measured hard colors. We 
have calculated the hardness ratios and photon fiuxes 
that we would expect to get from these energy bands for 
a variety of spectra and 0.5-8.0 keV luminosities using 
PIMMS and XSPEC. In Figure El we plot the colors and 
fluxes expected for power-law spectra with the dotted 
lines, and for a optically-thin thermal plasma with the 
solid lines. We have assumed a distance of 8 kpc and 6 x 
10^^ cm^^ of absorption from interstellar gas and dust. 
The median hard color for the Galactic center sources 
is —0.05. This corresponds to a r=1.5 power law or a 
kT— 1.7 keV plasma spectrum. 

The Galactic center sources in this survey are signifl- 
cantly softer than those from the deeper (limiting lumi- 
nosity of Lx = 2 X 10^^ erg s~^) catalog from the central 
20 pc of the Galaxy. The latter catalo g had a median 
hard color of —0.22 ijMuno et alJl200,^ . corresponding 
to a r=0 power law spectrum. This suggests that the 
more luminous X-rays sources are systematically softer. 
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Fig. 3. — The hard color plotted against the photon flux from each source. The symbol shapes indicate which observations the sources 
were identified in: blue squares for the Sgr B2 field, red stars fo r the Arches regio n, and green triangles for the extended shallow survey. 
For comparison, we also include the data from the central 25 pc IMuno et alJI20CT) using black circles. Open symbols indicate foreground 
sources, and filled circles those at or beyond the Galactic center (see text). Sources detected in only in the 3.3-4.7 keV band are assigned 
hard colors of —1; those only detected in the 4.7-8.0 keV band are assigned H R2=+l, and those detected in neither band are assigned 
H R2=—l.l. The uncertainties on the hard colors are significant. Sources with a probability of <50% of being detected have unreliable 
hard colors; those with a 50-90% chance of detection have aHR2 ~ 0.6; those with a 90-99% chance of detection have aHR2 ~ 0.4; and 
those with a >99% chance of detection have crHR2 ^ 0.3. Finally, we have plotted the colors expected for sources of varying luminosities 
at a distance of 8 kpc, and absorbed by 6 X lO^'^ cm"'^ of interstellar gas and dust. The dotted lines are for power-law spectra, and the 
solid lines for thermal plasma spectra. The sources with HR2>0.5 either have large uncertainties, large absorption columns, or both. For 
a fiducial conversion factor between photon flux and 0.5-8.0 keV luminosity, we assume a r=1.5 power law or a kT=7 keV plasma, and 
find that 10'^'' erg s~^ equals 6 X 10~^ ph cm^^ s^^. 



For a fiducial F = 1.5 spectrum absorbed by A'h=6 x 
10^^ cm"^), the photon fluxes can be converted to energy 
fluxes according to 1 ph cm~^ s~^ = 8 x 10~^ erg cm~^ 
s-i (0.5-8.0 keV). The de-absorbed 0.5-8.0 keV flux is 
approximately 3 times larger, so that for a distance D=8 
kpc, 10^'' erg s~^ equals 6 x 10^^ ph cm^^ s^^. For 
sources detected below 2.0 keV, we find that 1 ph cm^^ 
s-i = 2 X 10-9 erg cm'^ s'^ between 0.5-2.0 keV. The 
absorption for these sources is relatively small (< 10^^ 
cm~^), and therefore so is the correction to derive an 
intrinsic flux. 

3. RESULTS 

3.1. Spatial Distribution 

We examined the spatial distribution of the X-ray 
sources in order to determine how it compares to that 
of the ordinary stellar population observed in the in- 
frared. We were most interested in sources near the 
Galactic center, so we only considered those sources with 



soft colors HRO > —0.175. We also took care to select 
only sources that were bright enough to be detected over 
a large fraction of the survey. To do this, we derived 
maps of our sensitivity as described in Appendix A, and 
we examined only those sources that (1) were brighter 
than a well-defined flux limit, and (2) that were located 
at a position where the sensitivity was lower than that 
flux limit. Our flux limit was designed so that sources 
brighter than the limit have at least a 50% chance of 
being detected over a significant fraction of our survey 
(see Appendix A for a description of the Monte Carlo 
simulations that we used to calculate the detection prob- 
ability). We found that a flux limit of 3 x 10~^ ph cm~^ 
s~^ (equivalent to 5 x 10'^^ erg s~^ [0.5-8.0 keV] taking 
D=8 kpc, Nn = 6 X lO^^cm"^, and a r=1.5 power law) 
provided the largest num ber of sources, 321 out of the 
1899 sources in Table IbTI 

We then computed the surface density of sources as a 
function of offset from Sgr A* and of the absolute val- 
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Fig. 4. — The distributions of point sources as a function of 
angular offset from Sgr A* (top panel), and the absolute values of 
Galactic longitude {middle panel) and latitude {bottom panel). We 
only considered sources that were brighter than Fx = 3 X 10~® ph 
cm~^ s~^, that had a 50% chance of being detected, and that lay 
in regions where the 50% detection threshold was lower than the 
above flux limit. The contribution of each source to the distribution 
was weighted by the inverse of the probability of detecting it. The 
surface density of stellar mass is plotted with the solid line, which 
has been normalized through a chi-squared minimization to match 
the surface density of X-ray sources. In all cases, the normalization 
implied that there were 4 X 10~^ X-ray sources with Fx > 3 X lO^'' 



ph cm ^ s ^ for every 1 Mq of stars. 



ues of the Galactic latitudes and longitudes, and plot- 
ted them in Figure 2| To account for the variations 
in our sensitivity over the survey region, the contribu- 
tion of each source to the distribution was weighted by 
the inverse of the probability of detecting it (calculated 
as described in A ppendix A). We used the catalog of 
IMuno et all l|2003|) to fill in the inner 8' around Sgr A*. 
We used the identical flux cuts in the Sgr A* field as 
for the rest of the survey, which allowed us to include 
another 42 sources. 

To compare the distribution of X-ray sources with that 
of the ordinary stellar population, we used models for the 
Galactic bulge and the central 150 pc that we re derived 
from infrared maps bv lLaunhardt et al. (2002), and the 
exponential model of the Galactic disk in Kent, Dame, & 
Fazio (1991). These mass models are accurate to about 



50%. We modeled the central 150 pc with two compo- 
nents. We assumed that the central 15 pc is dominated 
by a spherical cluster with a mass density profile 

Pc 



P - 



1 + {r/rcY 



(3) 



For r<6 pc, we use Pc=3.3 x 10^ Mq pc"'^, rc=0.22 pc, 
and n=2. For 6 < r < 200 pc, n=3, rc remains the same, 
and Pc is adjusted so that the function is continuous at 
r = 6 pc. The total mass of the central cluster is 6 x 
10^ Mq. 

The rest of the central 150 pc is dominated by a disk- 
like distribution with a mass density 

P = P<ir-""exp(-|z|/zd). (4) 

For r<120 pc, we take pd=300 Mq pc~^, n=0.1, and 
Zd— 4:5 pc. For 120<r<220 pc, we take n=3.5, leave Zd 
the same, and adjust pd so the function is continuous at 
r=120 pc. For 220<r<2000 pc, we take n=10, and treat 
the other parameters the same as above. The total mass 
of this nuclear stellar disk is 1.4 x 10^ Mq. 

We model the Galactic bulge as a tri-axial ellipsoid of 
the form 



1/ci 



(5) 
(6) 

(7) 



The axis defining x,y,z is rotated 15° in east and 
1° north from our line-of-sight. The parameters are 
aa:=1100 pc, aj,=360 pc, az=220 pc, c^ = 1.6, and C||=3.2, 
and pbuigc=8 Mqpc"'^. The total mass of the bulge is 
taken to be 10^° Mq. 

Finally, we model the Galactic disk as a simple expo- 
nential, 

p = poexp(-r/rd)exp(-|z|/zd), (8) 

where = 2.7 kpc, Zd — 200 pc, and po = 5 Mq pc^'^, 
so that the total mass of the disk is 10"'^^ Mq. 

We integrated the model stellar density along our line 
of sight toward the Galactic center, using a lower limit 
of 6 kpc because we excluded foreground sources from 
the profiles in Figure 01 and an upper limit of 10 kpc 
because sources beyond this distance will be heavily ab- 
sorbed and difficult to detect. We then compared these 
model surface densities to the observed surface density 
of X-ray sources through a linear minimization. We 
did not attempt to correct the surface density of X-ray 
sources to account for fact that the absorption column 
varies as a function of longitude, latitude, and along our 
line of sight in the image, because the uncertainty intro- 
duced by our failure to do so is smaller than the ~50% 
uncertainty in the mass model. The best-fit stellar sur- 
face densities are indicated by solid lines in Figure ^ 
The overall match is good, with /v < 1. However, 
within 2' of Sgr A* (4.7 pc in projection) the number 
of X-ray sources is 2a larger than that expected from a 
simple scaling of the mass distribution inferred from the 
infrared (1.3±0.4 sources arcmin"^ versus the predicted 
0.5 sources arcmin"^). This may be further evidence 
that X-ray sources are more conce ntrated near Sgr A * 
than ordinary stars are (see also iMuno et al.| l2005a(l . 
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TABLE 2 

Parameters of the log N — log 5 Distribution 



Field 


■S'lim 


Num. 


Area 


a 


No 


Pks 




10-*^ ph cm-2 s-i 


Sources 


(arcmin^) 




(arcmin~^) 




Sgr B2 


0.7 


48 


101 


1.7±0.2 


0.04 


0.99 


Radio Arches 


2 


28 


97 


1.1±0.2 


0.17 


0.69 


Shallow Survey 


6 


142 


2833 


l.SitO.l 


0.11 


0.37 


Sgr A* 


0.8 


232 


152 


1.4±0.1 


0.24 


0.20 



Note. - — The normalization of the log N — log S distribution, A'^o is listed for a fiducial flux of 
3 X 10~^ ph cm~'^ s~^, to match the spatial distribution in Figure 4. Pks represents the probability 
under a Kolgoromov-Smirnov test of seeing the observed difference between the observed and model 
distribution assuming that they are identical, so that very small values would indicate a poorer 
match. 



1 .000 



^ 0.010 





Sgr A*: o=1.4+/-0.1 




^ Radio Arches; c<=1.2+/-0.2 ; 




WGL Survey: a=1.4+/-0.1 : 




Sgr BZ: c<=l,7+/-0.2 







10-6 10-5 10-4 

Photon Flux s-') 

Fig. 5. — The number-flux distribution of sources in the Arches 
Region (red lines), the Sgr B2 region (blue lines), and the shallow 
survey (blue lines). The solid histograms have been corrected for 
the detection probability, and are only plotted for sources brighter 
than the flux limits used in modeling the distributions. The 
solid lines are the best-fit model distributions {N{> S) oc 5""), 
with the slopes indicated in the upper-right of the figure. The 
dash ed lines represent th e ba ckground AGN contr ibutions taken 
from lBrandt et alJ 120011^ and lManners et alJ Imi^ . 



The normalization of the fits imply that above a limit 
of 5 X 10^^ erg there are (4±2) x 10"'' X-ray sources 
per solar mass of stars, where most of the uncertainty is 
in the mass models. We compare this to the expected 
density of X-ray sources in §4. 

3.2. Number-Flux Distribution 

Spatial variations in the underlying population of X- 
ray sources could be identified by examining the relative 
number of faint and bright X-ray sources. Therefore, we 
have computed the cumulative log A^— log 5* distributions 
in four regions: the Arches field, the Sgr B2 field, and the 
general shallow survey (excluding the deep pointings), 
and the Sgr A* field (see also^vmo ct al. 2003). How- 
ever, in this case, we required that each source had at 
least a 90% chance of being detected, which was stricter 
than when we studied the spatial distribution. We mod- 
eled the differential number-fiux distribution using the 
method described in Murdoch, Crawford, & Jauncey 
(1973), with slight modifications described in Appendix 
B to use Poisson statistics and knowledge of the aver- 
age background rate. As for the spatial distribution, we 



defined a flux limit (S'lim) for each region above which 
we could securely detect sources over the largest possi- 
ble area. For the logiV — \ogS distribution, we added 
the criterion that the fiux limit allow us to measure the 
fiuxes of sources at the ba level. Our choice of flux limits 
is described in more detail in Appendix B. The flux lim- 
its, area covered by each survey, and number of sources 
that were brighter than and located at points where 
sources with S—Sum could be detected securely are listed 
in Table El 

In Table [3 we also list the best- fit slopes a, the nor- 
malizations Nq at S'o=3 X 10~^ ph cm~^ (the fit 
had to be extrapolated for the shallow survey), and the 
probabilities Pks that the model and observed distribu- 
tions match according to the KS-test. The increase in 
the normalization of the distributions is consistent with 
the radial distribution of point sources in Figure 01 We 
find marginal, l.Scr evidence that the log — logS dis- 
tribution is flatter in the Radio Arches region than in 
the Sgr B2 field or the survey as a whole, which would 
imply that the former contains a larger proportion of 
high-luminosity sources. Unfortunately, our constraints 
on the log — log S distribution are a bit poor, because 
the number of sources that meet the criterion of having 
>5(T flux measurements (roughly, > 25 photons) is small 
(see Fig. EJ. The best-fit distributions are plotted with 
solid hues in Figure [3 

In Figure El we also plot the expected distributions 
of background active galactic nuclei (AGN). To convert 
between the photon fluxes we use here and the 2-8 keV 
energy fiux that is reported in papers on the Chandra 
and XMM-Newton deep fields, we assume that AGN are 
observed through an absorption column of 12 x 10^^ cm~^ 
and that their spectra are described by a r=1.5 power 
law (the final results are not sensitive to a choice of F 
between 1.2 and 1.8). We find that a source with a 0.5- 
8.0 keV photon flux of So = 3 x 10"^ ph cm'^ s'^ in 
our survey would have a 2-8 keV fiux of 6 x 10^^^ erg 
cm~^ if it were located at high Galactic latitude. 
Based on the number-flux distribution in the deep fields, 
expected density of background AGN is given by A^(> 
S) = NoiS/So)-'^, where at So = 3 x 10_f_ph cm_2 s-\ 
iVo=0.007 source arcmin"^ in Brand t et alJ i200lD . and 
iVo=0.01 source arcmin"^ in iManners^t'ldri|2003f) . This 
normalization is <10% of that of the shallow survey in 
Tabled (see also Fig. El. We conclude that <130 of the 
absorbed sources in our sample should be extra-galactic. 
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4. DISCUSSION 

We have presented a catalog of X-ray sources with 
fluxes between a few 10"^^ and 10'^^ erg s'-^. The majority 
of the region was covered with tiled pairs of 12 ks observa- 
tions, in which the 50% completeness limit is 1 x 10'^^ erg 
s~^. Two deeper exposures were more sensitive. A 50 ks 
toward the Radio Arches was complete to 4x 10"^^ erg s~^, 
and a 100 ks observations toward Sgr B2 was complete 
to 1 X 10^^ erg s~^. The sensitivity of the Radio Arches 
observation was a factor of two poorer than would be 
naively expected based on the exposure time, because the 
Galactic center produces strong diffuse emission against 
which point sources are more difficult to detect. 

Our survey encompasses a non-trivial fraction of 
the mass of sta,rs in the Galaxy. Integrating the 
(jLaunhardt et alJl2002j) models for the stellar distribu- 
tion over our survey area, we find that it encloses a stel- 
lar mass of ~10^ Mq, or 1% of the Galactic value. In 
Table O we list the various classes of object that could 
be detected as X-ray sources in our image if they were 
located near the Galactic center. We also include pre- 
dictions for the total number of each class of sources en- 
compassed by the survey region, based on the references 
in the table. The estimates are based on the assumption 
that the star form ation rate is ^1% of the Galactic value, 
or 0.01 Mq yr~i (jFiger et al.ll2004j) . However, if the re- 
cent rate of star formation in the central 150 pc of the 
Galaxy is ^--^10% of the total Galactic rate, as is suggested 
by indirec t measuremen ts of the Lyman-alpha fl ux in 
the region l)Cox fc Laure iis 1989; Fig er et aljfl999() . then 
WR/0 stars, luminous pulsars, and high-mass X-ray bi- 
naries each could be an order-of-magnitude more numer- 
ous. Below we describe the considerations that went into 
developing that table, and some implications that the ob- 
served population of sources have for understanding the 
evolution of the accreting binaries that make up the ma- 
jority of our sample. 

4.1. Comparison to the Local Galactic X-ray Population 

We start by comparing the amount of X-ray flux per 
unit stellar mass from the point sources with Lx>5 x 
10^^ erg in our survey to that in the local Galaxy 
as identified by 'Sazon ov et alJ l)2005fl . For a cumula- 
tive numbcr-flux distribution of the form A'^(> Lx) = 
No{Lx/ Lxfl)^" , the specific luminosity produced by 



sources with L> 



< Lx < Ly 



L 



X,tot 



uNqLx.o 



a - 



1 



Lx,n\\n 

Lx,o 



^X.max 

Lxfl 



(9) 

where iVo is the normalization at a luminosity of Lx,o in 
units of sources per solar mass. For the Galactic center, 
we have found that a = 1.5 and that Nq = (4± 2) x 10"'' 
sources M©"! at Lx,o = 5 x 10^2 erg s'^ (Tab. 13 Fig.^l. 
So, the the specific luminosity of X-ray point sources 
with luminosities between Lx,min = 5 x 10^^ erg s~^ and 
ix.max = lO^'' erg s~i in our survey is £x tot = (5 ± 2) x 
10^^ erg s"^ Mq^^- In the local Galaxv. FSazonov et al.l 
(j2005ll) find that a w 1.2 and iVo = (6± 2) x 10"'* sources 
Mq~^ at ix,o = 2 X 10^° erg s~^ (where for K in their 
Eq. 5, iVo = K / and we have estimated the uncertainty 
based on that of total specific luminosity for sources with 
Lx < 10^** erg s~^). So, the specific luminosity of sources 



TABLE 3 

X-Ray Sources in the Galactic Center 



Object 


log(-Lx) 


Number 


Number 


References 




log{erg s-i) 


in GC 


Detectable 




CVs 


29.5-33.5 


10^ 


10^ 


[1,2] 


WR/O Stars 


31-34 


10^ 


10 


[3,4,5] 


Pulsars 


29.3-35 


10^ 


10 


[6,7] 


LMXBs 


30-39 


10^ 


10 


[8,9,10] 


HMXBs 


31-38 


10^ 


50 


[11] 



References. — [1] 'Vorbunt ct al! !l997f); [2] !Sazonov ot 

t005): [3] iFiecrctal. 2004); [4] Pollock (1987); [5] Bcrriiofcr ct 
997): [ el.Bcckcr Asclicnbacli 12002.1: [71.Cordcs fc Lazio 11997): 
Wiinands et al.[ 12002}): 191 Iliu fc lJ i2005l) : 1101 (Kong et all ^2002, 
[11] Pf aTTeT^ l20a 

Note, — We list order-of-magnitude estimates of the total popula- 
tion of various X-ray sources in our field, along with the number that 
should be detected in our survey. 



with luminosities between Lx,min = 5 x 10"^^ erg s ^ and 
ix.max = 1034 erg s-i is locally ^x.tot = (1.0±0.3) x 10^^ 
erg s~^ Mq~i. Therefore, the specific luminosities of X- 
ray sources with 5 x lO'^^ < Lx < 10^^ erg s~^ in the 
local Galaxy and in the Galactic center are consistent 
within their uncertainties. 

To understand the population of X-ray sources at the 
Galactic center, it is notable that in the local Galactic 
neighborhood all the X-ray sources with lO'^^ < Lx < 
10'^'^ erg s~^ (2-10 keV) are cataclysmic variables (CVs) 
with magnetic white dwarfs and orbital periods of several 
hours (intermediate polars; Sazonov et al. 2005). There- 
fore, it is conceivable that most of the X-ray sources with 
Lx < 10"^^ erg s~ ^ in o ur survey are magne tic CVs (see 
also lMuno et alJl2004HlLaycock et al.ll2005j) . We would 
expect CVs to have the same spatial distribution as the 
old (> Gyr) population of stars, which dominate the in- 
frared light from the Galaxy, and indeed the distribution 
of X-ray sources is identical to th e inferred dis t ributio n 
of stars (Fig. EJ. As described in llVIuno et alJ l|2004b(l . 
intermediate polars have particularly hard, intrinsicly- 
absorbed spectra that are consistent with those of the 
sources with Lx ^ lO'^'^ erg s~^ in Figure|21 

Few CVs have Lx > lO'^^, and only one CV has been 
observed a t Lx ^. 5 x 10'^'^ erg s~^ (^-^ P^'" 
burst; e.g.. lKing. Ricketts. fc Warwick! 1 197^ . so we ex- 
pect brighter sources to be more luminous objects such 
as LMXBs, HMXBs, WR/O stars in coUiding-wind bi- 
naries, or pulsars (Table |2J|. The excess of bright X- 
ray sources observed in the Radio Arches field (Fig. O, 
in which the Arches and Quintuplet cluste rs are strik- 
ing evidence recent active star formation (i Figer et alJ 
1999), can be explained if some of the bright sources are 
HMXBs, WR/O stars, or young pulsars. Such sources 
have short lifetimes, and so should be concentrated near 
regions of active star formation. Such sources also have 
softer spectra than intermediate polars in the 2-8 keV 
band, and so could explain why the sources brighter than 
^-[^q33 gj.g g-i FigureOare systematically softer than 
the faint ones. 

4.2. Comparison to Theoretical Models 

Several binary population synthesis calculations have 
been carried out in order to interpret the population 
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of X-ray sources in t he Galactic center described by 
iWang et all pOOl and lMuno et"!!! (2003). These mod- 
els outline how the numbers of each class of X-ray 
sources constrain various combinations of parameters in 
the binary evolution models and assumptions about the 
physics of syst ems accretiiig at lo w rates. 



For instance. IP fahl et "all l)2002j) su ggested that s ever al 
hundred of the X-ray sources in the iWanE et al.' ("2005) 
survey could be neutron stars accreting from the winds of 
>3Mo binary companions (see also iBelczvnski fc Taaml 
I2OO4') . The total numbers of wind-accreting neutron 
stars, and the fractions of systems with companions more 
and less massive than SMq, varies by a factor of a few 
depending upon the the magnitudes of the kicks im- 
parted to the neutron stars at birth. These theoretical 
predictions have motivated searches for infrared coun- 
terpa rts to the X-ray sources in the Galactic center sur - 
veys llLavcock e t al.' "20051 iBandvopadhvav et al.ll2005l) . 
However. ILiu fc Li ( 200l^suggested that accretion could 
be inhibi ted by t he magnetospheres of the neutron stars 
(see also iDavies fc Pringle 1981) at the l ow mass-transfer 
rates considered by iPfahl et all l)2002j) . so that wind- 
accreting neutron stars may not be luminous enough to 
be detected in our Chandra survey. 

There also is theoretical disagreement as to whether 
CVs contribute significantl y to the popula tion of X- 
ray sources in our survey. ILiu fc Lil l)2005D . who base 
their calculations on the binary evolution code of Hur- 
ley, Tout, & Pols (2002), suggest that CVs are not lu- 
minous enough to be detected in large numbers from the 
Galactic center. However, this disagrees with similar cal- 
culati ons b y Ruiter et al. (2005), who use the StarTrack 
code l|Belczvnski et aL.2005') . and predict that significant 
numbers of luminous CVs should be detectable from the 
Galactic center. The main difference between the two 
codes is in their prescriptions for calculating the rate of 
mass transfer, which leads ILiu fc Lil l|2005ll to predict 
m ass transfer rates ^ 100 times lower than those used 
bv lRuiter et alJ l)2005|) for identical systems with orbital 
periods of several hours (A. Ruiter, private communica- 
tion) . Systems with orbi ts of severa l hours have the high- 
est accretion rates ( Patterson"1984l see also Howell, Nel- 
son, & Rappaport 2001) and therefore are the most lu- 
minous in X-rays, which makes them the most important 
contributors among CVs to our survey. Our comparison 
with the local Galactic population of X-ray sources sug- 
gests that CVs are indeed both numerous and luminous 



image (see also, e.E..lVerbunt et alJl997 


lEzuka & Ishidal 


ll999ilSuleimanov. R.evnivtsev. & Rittei 


I200fj). so our re- 



suits could be taken as further, indirect evidence i n sup- 
por t of the pres cription used bv 'Bel czvnski et aTl l)2005fl 
and iRuiter eTal. 12005) . 

Fina llv, IBelczvnski & TaamI l|2004() and iLui fc Lil 
l)2005D predict that there are a few thousands LMXBs 
consisting of a neutron star accreting from a white dwarf 
in our survey region. In these models, most of the neu- 
tron stars in these systems form through the accretion- 
induced collaps e of an ONe w hite dwarf. Under the 
assumptions of ILiu fc Lil l)2005tl . a few percent of the 
LMXBs are persistently bright enough to be detected 
in our survey. Moreover, most of these LMXBs should 
be transient, so if one assumes a standard duty cy- 



cle of ^^1%, there should be >50 LMXBs in outburst 
with Lx > 10'^^ erg s~^ in the field at any given 
time (and thousands in the Galaxy). In contrast, there 
are only two persistent LMXBs this lumi nous in the 
field, IE 1743.1-2843 and IE 1740.7-2942 i Wang et alJ 
'2002), and several decades of occasional X-ray observa- 
tions have only revealed about a dozen transients with 
Lx > lO'^^ erg s~ ^ (and on ly ^150 in the rest of the 
Galaxy; see, e.g.. iLhi. van Paradiis. fc van den Heuvel 
I2OOII: iWiinands et a l." 20051. The production of a large 
number of these transient LMXBs appears to be a com- 
mon feature of models in which neutron stars can form 
through accretion-induced collapse (e.g., §5.6 of Iben, 
Tutukov, fc Yunge l son 1 995). Under t he mo dels of 
IBelczvnski fc Ta am' (20Q3) and ILiu fc Lil |2005|) . possi- 
ble ways to accommodate the small number of bright 
transients in our field are to assume that the efficiency 
with which the envelope of a star can be ejected dur- 
ing the common envelope phase is low so that many of 
the binaries merge rather than forming LMXBs, or to 
assume that the accretion-induced collapse of an ONe 
white dwarf does not form a neutron star. 

5. THE FUTURE 

Further progress in understanding the natures of the 
X-ray sources near the Galactic center, and the con- 
sequential constraints on the parameters input into bi- 
nary evolution and population synthesis models, will 
be acquired through multi-wavelength observations of 
the region. X-ray observations will a dd to the popula- 
tion of transient LMXB s in the field (Wang et al. "2002 
i g^kano et al . 2005: Mu no et al.l Ii2005a: Wiinands et aD 
l2005i e.g.,). Comparing infrared and X-ray surveys 
will reveal individual e xamples o f WR/ binaries and 
HMXBs ( Bandvopadhv av et alJ2' 005: Mu no et al.l20fl5hl 
e.g.,). The first pulsars near the Galactic center may be 
found by comparing radio and X-ray surveys (see also 
iCordes fc Laziolll997^ . 

Finally, we note that this catalog will be improved 
upon greatly in the next year, because a series of deep, 
40 ks exposures of roughly half of the survey area have 
been approved for the 2006 Chandra observing cycle. 
We expect to increase the number of X-ray sources de- 
tected by a factor of ^5, and improve the uncertainties 
on the positions of the sources in the shallow survey from 
wl" to <0'.'5. In the meantime, the current catalog pro- 
vides the best available sample for studying the spatial 
and number-flux distributions of X-ray sources near the 
Galactic center, and for identifying their counterparts at 
other wavelengths. 
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Fig. A1. — The probability of detecting a source as a function of count rate (y-axis) and offset (line styles). The detection probability 
depends upon the survey region because of variations in the diffuse X-ray background; the longer observations are background-limited. 
The probability at a given count rate decreases as the offset increases, because the size of the PSF increases. We note that the curve for 
the 8—10' offset in the shallow (WGL02) survey was not used much, because sources in that survey almost always were found closer to the 
aimpoint of at least one of the tiled pointings. 

APPENDIX 
COMPLETENESS OF THE SURVEY 

In order to study the spatial and luminosity distributions of the X-ray sources, we need to calculate the limiting 
flux at which we can confidently d etect source s as a function of position in the survey. We performed synthetic star 
tests following the basic method of lBauer et al . (2004), with modifications to account for differences in our treatment 
of the photometry and the more complex layout of our survey. 

In order to produce ~10^ synthetic stars, we simulated 1000 exposures of ObsID 2287 (12 ks exposure), and 100 each 
for ObsIDs 945 (50 ks exposure), and 944 (100 ks exposure). For each observation, we removed events from within a 
circle circumscribing 92% of the energy of the PSF around each detected source. We then created images from the 
resulting event lists, and filled the "holes" in the image with a number of counts drawn from a Poisson distribution 
with a mean equal to that of a surrounding annulus. We drew fluxes for synthetic stars from a power-law distribution 
A'^(> S) cx S~°' with a slope a = 1.5, normalizations and flux limits chosen to match the numbers and intensities of 
sources in each region (e.g., Tab.l^J- We converted these fluxes to count rates using an exposure map generated for 
photons with E = 3 keV, and then drew net numbers of counts from Poisson distributions with those mean count 
rates. Each point source was assigned to a random position in the image. Next, we obtained a model image of the PSF 
created by mkpsf when computing the photometry for the the nearest detected source in the image, and used the PSF 
as the probability distribution to simulate the 2-dimensional image of the counts. This was placed into the composite 
image. Finally, we search the synthetic image for point sources using wavdetect, and tabulated which sources we 
placed into the image were detected by our search algorithm. 

We display the fraction of sources that are detected in 12, 50 and 100 ks as a function of input flux and offset from 
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Fig. A2. — The number of counts required to detect a source in 50% (dashed line) and 90% (solid line) of trials, as a function of offset 
from the center of the cluster. In order of decreasing number of counts required, the red line is for the Arches region, the blue line for Sgr 
B2, and the green line for the shallow survey. 



the aim point in Figure lAll The probabihty of detecting a source obviously increases with larger flux, because there 
is more signal-to-noise. For a given flux, the probability of detecting a source is also larger at small offsets from the 
aim point, because farther from the aim point counts are spread over a larger area by the PSF. This latter trend 
can also be seen in Figure IA2I where we have plotted the number of counts above which a source will be detected 
50% and 90% of the time as a function of offset from the aim point. As expected, longer observations enable us to 
reliably detect sources with lower fluxes. However, the total background counts from diffuse X-ray emission toward 
the Galactic center also increases for longer exposures, so more counts are required to reliably detect a source. As a 
result, a factor of 5 increase in exposure from 10 to 50 ks yields only a factor of 3 increase in sensitivity, which is by no 
means linear, but it is still a faster improvement than the t^^/^ trend that would be expected if we were background 
limited. 

We used the above simulations to make a map of the sensitivity of our survey. For each observation, we generated a 
sensiti vity map by (1) computing the offset of each pixel from the aim point, (2) computing the count threshold from 
Figure IX2l and (3) dividing the count threshold by the exposure map (in units of cm^ s) to obtain a flux. Then, to 
create a composite map for the full survey, for each pixel we recorded the lowest value of the flux threshold from all of 
the maps with exposure at that pixel. We note that using this method, localized enhancements in the diffuse emission 
that decrease our sensitivity have been averaged over offset from the aim point. Properly accounting for all of the 
variations in the diffuse emission would require that we carry out our Monte Carlo simulations for each exposure, which 
would be time consuming and would not change our results significantly. The resulting map is displayed in Figure IXSl 
We als o ha ve tabulated the area over which our survey is sensitive as a function of limiting flux, and displayed that in 
Figure El 

We have also used the simulations to compute the probability of detecting each source. For each source, we deter- 
mined the offset from the aim point of the most sensitive exposure, and determined the count threshold at that offset 
from Figure We then divided that count threshold by the exposure map from the most sensitive observation at 
the location of the source to obtain a flux threshold. Again, this technique ignores local background variations that 
introduce systematic errors in our computed thresholds for individual sources, but these errors should average out 
when considering large numbers of sources. The median probability of detecting a Galactic center source was 57%. 
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This is because the photons from faint sources are often lost to the wings of the PSF, and because the number-flux 
distribution is quite steep (see H3.2(l . so many faint sources are only detectable when Poisson fluctuations result in 
larger observed counts (this produces the Eddington bias when computing number-flux distribution). 

The sensitivity maps and detection probabilities were used in computing the spatial distribution in H3.ll and flux 
distributions in i j3.2l We also repeated the above process for the combined image created from 625 ks o f expo sure on 
the 17'xl7' field around Sgr A*, so that we could compare the results from the catalog in Muno et al.i l|200.lD t o th e 
current one. We find that a source can be detected confidently with the fewest counts in the shallow survey fFig. IA2|I . 
because the longer observations were background limited. Of co urse, the observations were still more sensitive when 
the detection threshold was considered as a function of flux fFig. IA4|l . 

MODELING THE NUMBER-FLUX DISTRIBUTION 

We modeled the un-binncd numbcr-flux distributions using the technique described in iMurdoch et all l)1973f) . We 
assumed that the cumulative number-flux distribution could be described as a power law A'^(> S) oc S~°', where S 
is the net number of counts for a source, and that the number of observed counts could be described by a Poisson 
distribution with mean rate S + B, where B is the average number of background counts in the source extraction 
region (e.g., i?=4.3 for the shallow survey). The likelihood of observing sources with a distribution of total counts Ci 
(which are all integers) is then: 

roo -{S+B)( c , f>\C, 

5]lnP(a)^^PdeMif-W ±i±J^iVo^-("+i)d5, (Bl) 

where K is normalization over the range of count rates under consideration (Ci, Cu), 

P,,,,£:^!l!^;^±^7Vo5-("+i)d5. (B2) 

The normalization of the power law Nq drops out of Eauation lBll and so it was derived by setting the normalization 
K for the best-fit a equal to the observed number of sources. Finally, we compared the observed cumulative flux 
distribution to the model distribution Pi using a KS-test, to establish whether our power-law model is consistent with 
the data. 

Caution needs to be used in exercising these equations, as is described in detail in IMurdoch et alJ ljl973|) . First, 
the integrals over S will diverge as S approaches unless a the co unt rate from a source is inconsistent with the 
background rate B at the «5cr level (IMurdoch et alJll973l:lWand2004|) . Therefore, we restricted our analysis to sources 
brighter than the 5cr detection threshold for the average background level in each region. The mean background and 
count thresholds are listed in Table El In doing so, we were able to ignore the negligible contributions to the above 
integrals from beyond S « {d) ± 12{C,y/^. Second, our source-detection algorithm was designed primarily to reject 
false positives, and we find that false negatives occur for «20% of sources with count rates at the 5(t level above 
background. We use the factor Pdct,i in the above equations to account for the probability of detecting a source with 
a count rate Ci, which we determined from our Monte Carlo simulations. Third, in order to avoid being biased by 
bright sources detected in regions with poor sensitivity, we only considered regions of the image in which there was a 
>90% chance of identifying a source with a flux equal to the 5a detection threshold. We list in Table |2 the number of 
sources that met these criteria and the area over which we were sensitive. Finally, we note that we used average values 
for the background and the detection probability, even though both varied signiflcantly over the regions covered by 
each set of observations. The total counts Ci used above were the sum of the net counts derived from our photometry 
and the average background, which was then rounded to the nearest integ er. This was necessary to ensure that the 
integrand in Equations Bl and B2 were monotonic functions of Ci (see also lWanell2004l) . 

Our approach to modeling the logA^ — logS* distribution takes into account the possible Eddington bias, although 
it ignores the vast majority of faint sources. Different approaches are possible, and should yield similar results. For 
instance, Bauer et al. ( 2004) used Monte Carlo simulations to estimate the corrections to the flux required to offset both 
the Eddington bias and biases introduced by their method for deriving the photometry for each source, assumed that 
the photometric uncertainties were negligi ble, and analytically c omputed a from a maximum-likelihood distribution 
using corrected fluxes and the equations in IMurdoch et alJ l)1973j) . We have not implemented this technique, because 
our procedure for computing the photometry for each source was computationally prohibitive to incorporate it into our 
Monte Carlo simulations. We also note that their resulting slope could be somewhat biased, because they have assumed 
a log N — log S distribution in computing the flux correction, thereby pre-determining the effect of the Eddington bias. 

iWand l)2004() has presented a method that is almost equivalent to ours, in which he applied a redistribution matrix 
to convert a model distribution into an observed distribution, and then foun d best-fit parameters for the model using 
the chi-squared and Cash minimization techniques implemented in XSPEC (lArnaud et al.llT'996D . If the model and 
un-binned fluxes are compared using the Cash statistic (Cash 1979), the techniques are equivalent, although the use of 
a response matrix converts the integral in Eauation lBll into a sum. Using binned data and a chi-squared test obviously 
requires enough sources per bin that their numbers are approximately distributed as a Gaussian. That technique has 
the advantage of using tools that X-ray astronomers are familiar with, although it is conceptually more complicated 
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than our method. 
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TABLE Bl 

Catalog of Point Sources within 2°x0.8° of the Galactic Center 



Source 


ra dec 


Unc. Obsid Offset Tcxp 


C'nct 


^dct 


HRO 


HR2 


Fx 10"^ Flags 


(CXO J) 


(J2000) 


(arcsec) (arcmin) (ks) 










(ph cm~^ ^~^) 



174204.8 ■ 
174206.8 ■ 

174208.5 ■ 

174210.4 ■ 

174216.0 ■ 

174216.1 ■ 
174217.8 ■ 
174218.6 
174219.2 
174220.1 

174220.1 
174221.3 ■ 
174221.3 
174222.8 ■ 

174223.3 ■ 

174224.2 ■ 

174224.6 ■ 

174227.5 ■ 

174228.3 ■ 

174228.4 ■ 

174228.5 - 
174230.0 
174230.3 

174230.6 ■ 

174231.4 ■ 



295003 
293634 
294436 
293639 
293756 

293732 
293715 
293931 
294333 
293526 

293905 
294250 
294647 
294118 
293950 

293412 
294333 
292602 
294157 
293431 

293736 
293949 
294713 
294320 
294336 



265.52017 
265.52844 
265.53574 
265.54358 
265.56689 

265.56723 
265.57425 
265.57761 
265.58029 
265.58389 

265.58414 
265.58887 
265.58902 
265.59518 
265.59735 

265.60095 
265.60275 
265.61487 
265.61816 
265.61850 

265.61883 
265.62524 
265.62646 
265.62769 
265.63089 



-29.83436 
-29.60947 
-29.74361 
-29.61107 
-29.63249 

-29.62580 
-29.62091 
-29.65884 
-29.72591 
-29.59068 

-29.65145 
-29.71409 
-29.77988 
-29.68858 
-29.66397 

-29.57001 
-29.72592 
-29.43403 
-29.69925 
-29.57547 

-29.62687 
-29.66380 
-29.78708 
-29.72239 
-29.72681 



1.5 
1.0 
1.0 
1.0 
0.9 

0.9 
0.9 
0.9 
0.9 
0.9 

0.9 
0.9 
1.1 
0.9 
0.9 

1.0 
0.9 
1.2 
0.9 
1.0 

0.9 
0.9 
1.1 
0.9 
0.9 



2289 
2290 
2290 
2290 
2290 

2290 
2290 
2290 
2290 
2290 

2289 
2290 
2286 
2290 
2290 

2290 
2283 
2283 
2290 
2290 

2290 
2290 
2283 
2290 
2290 



11.3 
6.4 
6.5 
5.7 
3.9 

4.2 
4.1 
2.7 
4.1 
5.3 

2.6 
3.3 
6.9 
2.0 
1.7 

6.2 
3.6 
8.1 
1.8 
5.7 

2.7 
0.4 
7.0 
3.1 
3.4 



11.6 
11.6 
11.6 
11.6 
11.6 

11.6 
11.6 
11.6 
11.6 
11.6 

11.6 
11.6 
11.6 
11.6 
11.6 

23.2 
11.6 
11.6 
23.2 
23.2 

23.2 
23.2 
11.6 
23.2 
23.2 



15.4 



+8.0 



4.2 
8.2 
14.4 
8.7 



6.9 

+4.1 
2.7 
+5.5 
4.0 
-6.8 
-5.7 
+5.0 
4.5 



11.6 
10.6 
6.7 
4.5 
4.5 



+5.8 
5.3 
+5.6 
5.0 
+4.3 
4.0 
+3.8 
3.0 
+3.9 
3.0 



io-9±5;3 



2.7^ 



7+2.7 
-2.3 

8.6+^-'-' 



3.7 
5.8 



4.8 
+3.2 
2.8 
3.9 
3.7 



12.0 
8.6 

66.5 
3.3 

9 



+7.1 
5.3 
+5.1 
4.4 
+ 13.9 
13.1 
+4.2 
3.0 
+5.7 
-5.4 



+5.8 
-5.7 
+5.1 
4.1 
+5.4 
4.7 
+9.2 
7.7 

< 8.7 



7, 

4.3 
8.6 
20.2 



0.00 


'^•''^-0.42 


-9.00 


60.2 


0.25 


-9.00 


1.00-0.76 


22.8 


0.27 


-1.00+050 


-9.00 


24.8 


0.95 


-9.00 


20+0.46 
U-^y-0.45 


64.7 


0.81 


n 15+0'™ 

U.10_Q 73 


-9.00 


25.9 


0.99 


-9.00 


R7+0.33 

-0.34, 


57.0 


0.96 


1.00-1.09 


_0 0^+0.58 


36.0 


0.94 


1.00-1.09 


n 20+0-80 


25.2 


0.17 


-9.00 


-LOO+O-''" 


14.1 


0.06 


1.00-0.89 


-9.00 


13.5 


0.99 


-1.00+0.55 


1-00-0.87 


36.9 


0.02 


-9.00 


-1.00+1-08 


7.7 


0.63 


1.00-1.13 


46+0.55 
"■^"-0.64, 


40.4 


0.34 
0.69 


-9.00 
-9.00 


-0 00+0-99 
-n 'ici+0.77 


14.2 
20.6 


0.53 
0.73 
0.89 


1.00-0.85 

-0 01+0 09 
'^•'-'^-0.97 


-0 02+0-59 
'-'•'^^-0.56 
-1.00+0-64 

-0 20+0-''6 
'-'•^"-0.65 


29.7 
22.3 
204.3 


0.11 


-9.00 


-0 12+1-" 
^^•^^-0.88 


10.1 


0.09 


-1.00+042 


1.00-1.33 


15.4 


0.75 


-9.00 


1.00-1.13 


20.3 


0.10 


-9.00 


-9.00 


9.3 


0.17 


1.00-0.92 


-1.00+0-51 


24.3 


0.99 
0.00 


07+0 62 
'^■'^'-0.54 


-1.00+0.98 


46.9 
8.9 



sf 



sf 

f 



Note. — The full table -will be available online, and is and is on http : //astro .ucla. edu/inmuno/ sgra/ shallow_survey_catalog.txt 



